The Differential Equation Method: 
evaluation of complicated Feynman diagrams 

A.V.Kotikov 

Particle Physics Laboratory, Joint Institute for Nuclear Research 
141980 Dubna, Russia. 

Abstract 

We discuss a progress in calculation of Feynman integrals which has been done with help of the 
Differential Equation Method and demonstrate the results for a class of two-point two-loop diagrams. 



The idea of the Differential Equation Method (see fljl-lll) ( see a reviews in H ) : to apply the integration 
by parts procedure |^ to an internal n-point subgraph of a complicated Feynman diagram and later to 
represent new complicated diagrams, obtained here, as derivatives in respect of corresponding masses of 
the initial diagram. 

The integration by parts procedure ||], for a general n-point (sub)graph with masses of its lines 

1711,1712, ...^mn, line momenta pi , p2 =Pi—pi2,Pn = Pi — Pin and indices ji,j2, ■■,jn, respectively, has 
the following form: 
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where cu = pi 



are the propagators of n-point (sub)graph. 



Because the diagram with the index {ji -I- 1) of the propagator Ci may be represented as the derivative 
(on the mass m^), Eq.(^ leads to the differential equations (in principle, to partial differential equations) 
for the initial diagram (having the index ji, respectively). This approach which is based on the Eq.(n|) 
and allows to construct the (differential) relations between diagrams has been named as Differential 
Equations Method (DEM). For most interested cases (where the number of the masses is limited) these 
partial differential equations may be represented through original differential equation^ which is usually 
simpler to analyze. 

Thus, we have got the differential equations for the initial diagram. The inhomogeneous term contains 
only more simpler diagrams. These simpler diagrams have more trivial topological structure and/or less 
number of loops and/or ends H, ^. 

Applying the procedure several times, we will able to represent complicated Feynman integrals (FI) 
and their derivatives (in respect of internal masses) through a set of quite simple well-known diagrams. 
Then, the results for the complicated FI can be obtained by integration several times of the known results 
for corresponding simple diagrams |^ 

Sometimes it is useful (see [||) to use external momenta (or some their functions) but not masses as 
parameters of integration. 



The recent progress in calculation of Feynman integrals with help of the DEM. 

1. The articles Q and 



^The example of the direct application of the partial differential equation may be found in g). 

^In calculations of real processes (essentially in the framework of Standard Model) it is useful to use the relation (1) 
(at least, at first steps of calculations) to decrease the number of contributed diagrams (see ||^-[^ and |Q and references 
therein) . 
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Figure 1: Two- loop selfenergy diagrams. Solid lines denote propagators with the mass m; dashed lines 
denote massless propagators. 



a) The set of two-point two-loop FI with one- and two-mass thresholds has been evaluated by DEM 
(see Fig.l). The results are given on pages 2 and 3 and of some of them have been known before (see 

). The check of the results has been done by Veretin programs (see discussions in ^ and references 
therein) . 

b) The set of three-point two-loop FI with one- and two-mass thresholds has been evaluated (the 
results of some of them has been known before (see ^)) by a combination of DEM and Veretin programs 
for calculation of first terms of FI small-moment expansion (see discussions in |ll[ and references 
therein) . 

2. The article ||l|: 

The full set of two-point two-loop onshell master diagrams has been evaluated by DEM. The check of 
the results has been done by Kalniykov programs (see page 5 and discussions in JT^ , [l^ and references 
therein) . 

3. The articles 

The set of three-point and four-point two-loop massless FI has been evaluated. 



Here we demonstrate the results of FI are displayed on Fig.l. 

We introduce the notation for polylogarithmic functions ||l5| : 

TW^ c (-1)°+" /■ Mog°(t)log''(l-zt) 
Llo(z) = Sa-l,l[Z), ba+l,b[z) = ^, ^, J dt. 
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We introduce also the following two variables 



q2 _ i-yi7(i^ 

1 + V^/(z-4) 

Theng 

• /i - log2(-z) log(l - z) - 21og(-z)Li2(z) + 3Li3(z) - 65i,2(z) - log(l - z) (^(2 + 2Li2(z) ) , 

q^-h = 2C2 log(l + z) + 2 log(-z)Li2(-z) + log2(-z) log(l + z) + 41og(l + z)Li2(z) 
- 2Li3(-z) - 2Li3(z) + 251^2(2') - 451^2(2) - 45i,2(-z) , 

■ I12 = Li3(z) - 6C3 - C2 logy - ^ log^ y - 41ogy Li2(y) + 4Li3(y) - 3Li3(-y) + ^Usi^y^) , 

■ /i3 = -6Si,2{z) - 21og(l - z) (^C2 + Li2(z 

g2 . /i4 = log(2 - z)(^log2(l - z) - 21og(-z)log(l - z) - 2Li2(z)j - ^log^(l - z) - 2C2log(l - z) 
+ log(-z) log2(l -z)- ^1.2(1/(1 - z)2)+25i,2(l/(l - z))+25i,2( - 1/(1 - + i log^ y 
+ log2 y (^2 log(l + y2) _ 3 log(l -y + y^)^- 6C3 - Li3(-y2) + ^Li3(-y3) - 6Li3(-y) 

+ 2 logy (^Li2(-y2) - Li2(-y3) + 3Li2(-y)^ , 

q^ ■ /i5 - 2Li3(z) - log(-z) Li2(z) + (2 log(l - z) - ^ log' y (^8 log(l - y) - 3 log(l -y + y')^ - 6(3 
+ I log' y - ^Li3(-y=^) + 3Li3(-y) + 8Li3(y) + logy (^Li2(-y3) - 3Li2(-y) - 8Li2(y)^ , 

9' • /i23 = -C2 (^log(l - z) + logy^ - 6C3 - ^ log(l - y + y2) fog^ y + Li3(-y3) - 9Li3(-y) 

-21ogy (^Li2(V)-3Li2(-y)) , 
g' • /i25 = -2 log' y log(l - y) - 6C3 + 6Li3(y) - 6 logy Li2(y) , 

«' • /1234 = -6C3 - 12Li3(y) - 24Li3(-y) + 8 logy (^Li2(y) + 2Li2(-y)^ + 21og2 y (^log(l - y) 

+ 21og(l + y)V (2) 



Here we demonstrate the results of FI are displayed on Fig. 2. 

We consider here the following master-integrals in Euclidean space-time with dimension _D = 4 — 2e: 



^We would like to note that the coefficients of expansions of the results in respect of 2 are very similar (see ]16| ) to 
results for the moments of structure functions of deep inelastic scattering. 
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ONS{IJ}{i,j,m) EE R-^ I d^kP^'\k,Im) P'^^\k-p,Jm) 
J{JJIC}{i,j,k,m) EE K 
Y{JJJCC}{i,j,k,l,m) 



¥{ABIJK.}{a,b,i,j,k,m) 



where 



d'^kid^k2P^\ki,Im)P^iHki - k2,Jm) P^^\k2 - p^lCm) 
R-'^ j d^kid^k2P^''>{k2-p,Im) 
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the normahzation factor l/(27r)^ for each loop is assumed, and A,B,T,J^,IC ~ 0, 1. 

The finite part of most of the F-type master- integrals can be obtained from results of Ref . ||] in the 
limit z ^ 1. FlOlOl and Flllll have been calculated in Refs.|l^, |l^, respectively. Instead of the 
usually taken FOllOl integral |l^ we consider Jill as master integral. We recall the results of all 
master integrals for completeness. The last master integral FOOlll has been found in 

The finite part of the integrals of V- and J-type can be found in Refs.pot. The calculation of the e e 
(e^) parts of master integrals of this type have been performed by DEM. 

The results for F-type master-integrals are follows: 



F{ABIJIC}{1, 1, 1, 1, 1, m) ^ aiC(3) + aa^Sz + 03^^(2) + C'(e), 

V3 



(3) 



and the coefficients {oi} are given in Table I: 
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where ^ 



52 = ^Cl2 {-\ = 0.260434137632 ■ 
9\/3 V3/ 



Here we used the rr? — ie prescription. The results for the remaining master integrals are the following 
ones: 

vd^r^Ki, 1,1,1.,,.) . ^4(5-^),^4c,.)-4^-fs. + ^,„3 

+ y + ^2C(2) - 63C(3) - 12— - 63^2 + b^a2) ln3 + -^hS2^ 



63 

'a 



^52ln3 + 4-^ln3- ^-^In^S- ^^C(2) 



^^ln3-i^, 
x/3 2 V3 



2lLs3(2f) 



2 v^^' ' 2 V3 



where the coefficients {bi} are listen in Table II 0: 

''The results for the master integral VlOOl had a little error (see l22j) 
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Jlll(l,l,l,m) = 
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52C(2) + 48C(2)ln2-28C(3) 



0{e 



2(l-2e)(l-3e)U2 3^^^^ 



'7! 
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where [ffSl 
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^6* and Ls.ci — ^ -2.14476721256949- 



(4) 



(5) 



81 
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(6) 



(7) 



The above resuhs were checked numerically. Pade approximants were calculated from the small 
momentum Taylor expansion of the diagrams . The Taylor coefficients were obtained by means of the 
package |2^ with the master integrals taken from ||2^ . Further we made use of the idea of Broadhurst 
p6[ to apply the FORTRAN program PSLQ ||2^ to express the obtained numerical values in terms of a 
'basis' of irrational numbers, which were predicted by DEM. 

Let us point out that the numbers we obtain are related to polylogarithms at the sixth root of unityQ 
and hence arc in the same class of transcendentals obtained by Broadhurst in his investigation of 
three-loop diagrams which correspond to a closure of the two-loop topologies considered here. 
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^For the results obtained in expansion, however, the arguments of polylogarithms have other values (see Esj). 
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Figure 2: The full set of two- loop self-energies diagrams with one mass. Bold and thin lines correspond 
to the mass and massless propagators, respectively. 
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